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Abstract 



We study a statistical model defined by a conformally invariant distribution of 
overlapping spheres in arbitrary dimension d. The model arises as the asymptotic 
distribution of cosmic bubbles in d + 1 dimensional de Sitter space, and also as the 
asymptotic distribution of bubble collisions with the domain wall of a fiducial "obser- 
vation bubble" in d + 2 dimensional de Sitter space. In this note we calculate the 2-, 3-, 
and 4-point correlation functions of exponentials of the "bubble number operator" an- 
alytically in d = 2. We find that these correlators are free of infrared divergences, 
covariant under the global conformal group, charge conserving, and transform with 
positive conformal dimensions that are related in a novel way to the charge. Although 
by themselves these operators probably do not define a full-fledged conformal field 
theory, one can use the partition function on a sphere to compute an approximate 
central charge in the 2D case. The theory in any dimension has a noninteracting limit 
when the nucleation rate of the bubbles in the bulk is very large. The theory in two 
dimensions is related to some models of continuum percolation, but it is conformal for 
all values of the tunneling rate. 
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1 Introduction 

De Sitter space plays a central role in cosmology. In the standard model, the first moment 
after the big bang was a period of inflation in which the universe expanded many times over 
during an approximately de Sitter phase. Recent observations indicate that the expansion 
is currently accelerating, indicating the presence of a dark matter component which is most 
economically explained by a positive cosmological constant — in which case the future of our 
universe is a de Sitter phase. Moreover, if string theory is the correct description of nature, 
we may expect that our entire observable universe is a bubble expanding in an eternally 
inflating false vacuum de Sitter space [HE]- In such a model bubbles of the true vacuum 
nucleate in the false vacuum, and their walls accelerate outwards under the gravitational 
pressure due to the differing values of the vacuum energy inside and out. Observers may 
form inside these cosmic bubbles, which viewed from the inside appear to be negatively 
curved expanding Robert son- Walker cosmologies |3J. The cosmology inside will be affected 
in interesting ways both at its "Big Bang" jl] and later, by collisions with other bubbles (see 
e.g [5]). If this model is correct, understanding de Sitter space becomes even more crucial. 

But de Sitter space has proved maddeningly difficult to understand. In eternal de Sitter, 
points at fixed comoving distance always become causally disconnected at late times. As a 
result, correlation functions of points separated by more than a single Hubble length are not 
observable, and so one cannot define observables (like an S-matrix or boundary correlation 
functions in anti-de Sitter space) using well-separated points in space [6] . Moreover, because 
of the thermal nature of de Sitter and its finite entropy one expects correlators of operators 
inserted at timelike separated points to be quasi-periodic functions of the time separation, 
which neither converge nor settle into any predictable pattern (and in fact eventually produce 
fluctuations to every state consistent with the conservation laws) [7J |S]. Attempts to define 
a dS/CFT correspondence indicate that if such a dual theory exists, it cannot be unitary 
[10] . In the larger "multiverse" of the string theory landscape, one would like to understand 
how to average over distributions of cosmic bubbles so as to compute the expected values 
of cosmological observables visible to those living inside them. Much interesting work has 
been done on this problem; see e.g [UJ for a review. However, this too has been plagued by 
infinities and the non-uniqueness in the choice of measure. 

Here, we analyze the correlation functions of a putative dual CFT for eternal inflation. 
Our philosophy is that if one can find a consistent and well-defined set of quantities in 
eternally inflating spaces, this may lead to a solution of many of the problems above. To 
begin, consider an eternally inflating spacetime in d+1 spacetime dimensions in which at least 
one interacting scalar field is undergoing tunneling and forming bubbles. If one takes a d- 
dimensional constant time slice across the spacetime at late time, the slice will contain many 
casually disconnected regions and many vacuum bubbles, as shown in figures [T] and |2j Each 
bubble appears initially with small size, but as time passes its radius grows, asymptoting 
to a finite comoving radius determined by the conformal time at which it appeared (late 
appearing bubbles are smaller). 

We make some simplifying assumptions which allow us to compute the statistical distribu- 
tion of these bubbles on the <i-slice. (For a more general discussion, see [15].) This approach 
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Figure 1: A Carter-Penrose conformal diagram for bubbling de Sitter space. One can slice 
with a constant global time surface, which from d = 2 + 1 de Sitter would produce a 2-sphere 
tiled with bubbles (see Fig. [2]). Alternatively one can consider the sky as observed at late 
times by a "census taker" living in a A = bubble. Starting from d = 3 + 1 de Sitter, this 
will produce a very similar 2-sphere. 

has the advantage of simplicity, but the disadvantage that these bubble distributions are not 
observable. 

Any given observer can observe only a subset of the eternally inflating spacetime. Con- 
sider an observer inside a bubble of some type, the "observation bubble." The observation 
bubble collides with other bubbles that form through quantum nucleation in the false vacuum 
nearby, as shown in figure [TJ Each collision bubble will affect the observation bubble's wall 
inside a ball (a disk in the case of an observer 3+1 dimensional de Sitter). As time passes 
for an observer inside the observation bubble the collision appears as a point and grows at 
a rate which asymptotically approaches that of the observation bubble's wall itself. As a 
result the angular radius of the disk asymptotes at late time to some finite size if), which is 
determined in a simple way by the time of the collision (late collisions make smaller disks). 
If the observer has access to large conformal times — which requires that the cosmological 
constant inside the observation bubble be small — she can observe this asymptotic distribu- 
tion by, for example, its effects on her cosmic microwave background sky. Such an observer 
has been referred to as a "census taker" [T2] . 

We will investigate the statistics of the distribution of these bubbles and their collisions, 
focusing primarily on the case of 2D distributions. The distribution can be thought of 
as describing either the set of bubbles on a late time slice in 2+1 dimensional de Sitter 
space, or the set of collisions with an observation bubble in 3+1 de Sitter. As we will see, 
using the bubble distribution on these surfaces one can define a set of correlation functions 
which behave like primary operators of a 2D conformal field theory. The operators are the 
exponentials of a discrete number operator that counts the number of disks that overlap the 
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Figure 2: Artist's rendition of a global slice of bubbling de Sitter space at late times, or a 
census taker's sky. 

point where it is inserted. 

Using the distribution we compute the 2— ,3—, and 4-point functions of these exponential 
operators analytically and exactly. These turn out to be consistent with the hypothesis 
that these exponential operators are primary fields of a conformal field theory, except that 
the 4-point function exhibits a non-analyticity as a function of the location of the operator 
insertions. The theory is conformally invariant for arbitrary values of the dimensionless 
tunneling rate 7. The relation between the charge (3 of the exponential of the number 
operator e t/3N ^ and its scaling dimension is interesting and novel: 



We compute the central charge of this putative conformal field theory by evaluating its 
partition function on a sphere of radius R, and find a result proportional to the continuous 
parameter 7. 

We also discuss higher dimensional versions of the theory. The simplest such extension 
is a 3D version obtained from the statistics of bubbles on 3D global slices of 3 + 1 de Sitter. 
This 3D theory has a dimension zero number operator, and its exponentials again behave 
like primary fields with positive dimension. In fact, we can show that in any dimension 
d, the theory always contains a dimension zero number operator and exponential operators 
constructed from it with arbitrarily small positive weight. Since no such operator can exist 
in a unitary field theory in more than d = 2, these higher dimensional theories cannot be 
unitary. 

Finally, we will demonstrate that in the limit that the decay rate 7 — > 00, the theory in 
any number of dimensions becomes free: correlators of the exponential operators factorize 
onto products of 2-point functions in precisely the same way as vertex operators for a free, 
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Figure 3: A potential for a scalar field coupled to gravity that could produce bubble distri- 
butions of the type we consider. 



massless scalar in two dimensions. 



2 Conformal invariance 

The simplest case to consider is one in which the cosmological constant inside the bubbles is 
the same as that of the "false vacuum" outside, and where bubbles can nucleate both inside 
and outside other bubbles — always with the same decay rate 7, where 7 is defined as the 
dimensionless decay rate per unit Hubble time per unit Hubble volume. 

In order for the tunneling rate to be constant, all of the vacua should be identical. The 
simplest model is the potential drawn in figure |3j where the potential has a discrete shift 
symmetry relating the vacua to each other. Even once the potential is symmetric, there 
are in general nontrivial interactions between bubbles, such as correlations between their 
nucleation points and collisions between bubbles. We work in the noninteracting limit where 
the nucleation rate is constant, independent of the presence of other bubbles. 

With these assumptions one can compute the bubble distribution on a global time slice. 
The spacetime is de Sitter with metric 

ds 2 = {-dr? + dfijj) . (2.1) 

sin 77 

where dVt 2 d is the metric on a <i-sphere. The number of bubbles that nucleate in a conformal 
time interval drj is proportional to the spacetime volume available, 

dN = ^dVa+t = 7 — **-dfl d (2.2) 
sin 77 

where dV^ is the d— dimensional spacetime volume element and the dfl^ factor refers to the 
location of the nucleation point on the spatial slice. 

We want to characterize the distribution of bubbles on future infinity of de Sitter space, 
which is given by rj = in these coordinates. A given bubble nucleation nucleation will affect 
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Figure 4: A bubble nucleated at conformal time r\ covers an angular size ip = \r)\ on future 
infinity. 

a ball on future infinity. The domain wall of a bubble asymptotically approaches the future 
lightcone of the nucleation point. Light rays satisfy 

dip = drj (2.3) 

where ip is the angular radius of the lightcone. Therefore, a bubble nucleated at time rj has 
angular size ip = \r]\ on the conformal boundary, as shown in figure [ij The distribution of 
bubbles on the boundary is then 



dN = ~ . dili ■ (2.4) 



dip 
sin" 1-1 ip 

Previous attempts to define a theory on such slices can be found in [9j [13] . 

As mentioned in the introduction, we can also consider the distribution of bubbles which 



collide with a given "observation bubble." In [T5] we computed this distribution: the dis- 
tribution of collision bubbles on the boundary of the future lightcone of a point in 3+1 de 
Sitter space (i.e., the observation bubble's wall) after infinite time. The number distribution 
takes the following simple form: 



dN(ip, 9, 0) = —j-j sin 6d0d(f>, (2.5) 



47 dip 
3 sin 3 ip 

where ip G (0, it) is again the angular radius of the disk, and 9 and <p are the coordinates of 
its center on the boundary sphere (we have set the radius of the sphere to 1 for convenience). 



Notice that this distribution is identical in form to Eq. 2.4 with d = 2. An attempt to define 
a conformal field theory on this sphere can be found in [T4]. More generally, the number 
distribution of collisions on the boundary of an "observation bubble" in d + 2 dimensional 
de Sitter space is 

Q d+2 dip 

dN = 7^ . d+ i , dn d (2.6) 

ild+i sin ^ ip 
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Figure 5: Numerical simulations of the model in d = 2. Only one disk type is shown. Left 
pane: 7 = .1,5 /R = .01. Right pane: 7 = .5,5/R = .01. 

where Q n is the surface area of a unit n-sphere. 

This distribution turns out to have the remarkable property that it is invariant under 
global conformal transformations SO(d + 1, 1), which are Mobius transformations in d = 2. 
The easiest way to see this is to stereographically project the distribution to the plane. This 
projection maps spheres to spheres. In terms of the coordinates x l of the center of the sphere 
and the sphere radius r, a little algebra shows that the distribution becomes simply 

dN = ^d^ (2.7) 

where x % are the Cartesian coordinates on the plane. Note that this is also the small an- 



gular radius, small area approximation to the sphere distribution (2.4). Global conformal 
transformations on the sphere are generated by rotations plus special conformal transforma- 
tions. In the stereographic plane, special conformal transformations around the origin take 
a particularly simple form — they are simply the scalings x l — > Xx l and r — > Ar. Therefore 



the distribution (2.4) is conformally invariant, since it is manifestly rotation invariant, and 
one can always choose the origin of the stereographic plane to coincide with one of the fixed 
points of the special conformal transformation. As noted in |15j . this symmetry group is the 
set of Lorentz transformations on the point of nucleation of the observation bubble. 



2.1 Fractals and percolation 



2.7 in [To], where he commented that 



Benoit Mandelbrot considered the distribution Eq. 
in d — 2 it approximates the distribution of craters on the moon. One can easily calculate 
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the Minkowski or "box-counting" fractal dimension of various sets of points defined by this 
distribution. For example, the set of points not inside any bubble is a set of measure zero 
(because the set of spheres with radii in any given logarithmic interval covers a finite fraction 
of the volume). To compute the fractal dimension of this set, one chooses boxes of linear 
dimension e, finds the minimum number iV&(e) of boxes necessary to fully cover the set, and 
then the dimension is defined as cLf = lim e ^o( — din Nb/dlne). 

To compute this, note that the volume V(e) which remains uncovered by spheres of radius 
greater than e is the volume V(e + de) minus the volume covered by spheres of size between 
e and e + de: 

dV(e) = -C d e d dN (2.8) 
where Cd is the volume of a unit <i-sphere. The infinitesimal number of spheres is given by 



so we have the differential equation 



dN = 7^rV( £ ) , (2.9) 



% = -70,f (2.10) 



Integrating this equation gives V(e) = V e lCd . The number of cubical boxes of size e required 
to cover this area is simply N b (e) = V(e)e~ d , and therefore 



d F = d- 1 C d (2.11) 

(see [20J for this calculation in the context of bubbles in 3 + 1 de Sitter). Similarly one could 
compute the dimension of other sets, such as the set of points covered by exactly, or at most, 
k spheres. These will be scale invariant fractal sets of measure zero as well. The existence 
of these fractals is not surprising given the statistical scale invariance of the distribution. 

One may expect that these sets will undergo percolation transitions at special values of 
7 [16]. For example in D = 2, as one increases 7 from zero there should be a percolation 
critical point where the set transitions from connected along filaments to a disconnected 
dust. This will occur at some dp ~ 1. In d — 3 one expects two such transitions: from 
"ramified veils" to filaments at dp ~ 2, and from filaments to dust at some other dp ~ 1 



If in fact one can regard this as a type of percolation model, it has some novel properties. 
The most striking is that as we will see, the theory seems to be conformally invariant for 
all values of 7 rather than just at the critical points — or at least one can define operators 
that transform covariantly for all 7. Additionally, at least some correlation functions can be 
computed analytically and exactly up to at least the 4-point function. It would be interesting 
to see if quantities of primary interest for percolation (for example crossing probabilities) 
could be extracted using the techniques developed here, but we leave this question for future 
work. 

Models of a somewhat similar type have been considered in the past under the name 
"continuum fractal percolation," where "continuum" refers to the lack of an underlying 
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lattice, and "fractal" to some self-similarity in the distribution of disks (or other shapes — see 
e.g. [17J). Another model with some similarities is "Mandelbrot percolation," in which a 
square is subdivided into iV 2 smaller squares for some integer N, each of which is colored 
with a probability p, and then the remaining uncolored squares are subdivided and the 
process repeats. Interestingly, while this model does have a percolation transition, it is first 
order [18] — and hence not conformal even at the transition! Percolation in power-law disk 
distributions was considered in the context of networks in [T9] . 



3 Correlation functions 

In this note we will concentrate on the correlation functions of a field with multiple vacua, as 
pictured in figure |3j We label the minima by N, where N takes integer values. As mentioned 
earlier, we assume that the vacua are identical, so there is a shift symmetry iV —> N + 1. 
Starting from a given vacuum, the field can tunnel to the left or to the right; we call these 
events instantons and anti-instantons. 

When the forward lightcones of two spacelike separated nucleation points overlap, addi- 
tional physics is needed to determine the field profile in the overlap region. If the critical 
bubble size is small, one can think of these overlap regions as collisions between bubbles. It 
is then nontrivial and model dependent to solve for what happens in the future of a collision. 
In our model where the vacua are degenerate, the critical bubble is horizon size, so instead 
of collisions it is more accurate to think of bubbles nucleating on top of existing domain 
walls. We make the simplest possible assumption about the overlap regions: we assume the 
instantons satisfy superposition. In a region to the future of the nucleation points of iV + 
instantons and N_ anti-instantons, we assume the field is in the minimum N = N + — iV_. 
These simple assumptions allow us to calculate correlation functions explicitly, but it would 
clearly be interesting to perturb away from them by allowing interactions between nucle- 
ation points and nontrivial dynamics in the overlap regions. For small tunneling rate 7, the 
instantons are very dilute and one may expect that the interactions are unimportant. 

It is convenient to construct a partition function for the bubble distribution with which 
to compute expectation values. For convenience we will work on the plane, using the distri- 



bution (2.7). The partition function is 

Z = E £ A (P + + P-T J^J *** = ex P (7(P + + P-) jf % I d 2 ?) ■ (3.1) 



Each term in the sum corresponds to a configuration of n disks, with the kth disk centered 
at the point %k and with radius r^. The factor (p + + pJ) denotes the probabilities of the 
two possible types of nucleations, left-moving and right-moving. Because all of the vacua 
are identical by assumption, detailed balance demands p + — p_ = 1/2. To avoid infinities 
the integral must be cut off in both small and large disk sizes, although as we will see, well- 
defined correlation functions on the plane do not depend on the IR cutoff R. The factor of 
7™ is the appropriate weight for a configuration of n disks, given that 7 ~ e~ Sinst and that 
instanton interactions can be neglected. The partition function factorizes into instanton and 
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anti-instanton pieces, 

Z = Z + Z_ (3.2 

with 

" R dr 



Z+ = exp (j + ^ — J d 2 zj . (3.3) 

Symmetry determines 7+ = 7- =7/2. 

Because the distribution is invariant under Mobius transformations, one expects corre- 



lation functions of well-defined operators computed using the partition function |3.1| to be 
Mobius covariant. The subtlety arises from the cutoffs — divergent correlators will not trans- 
form simply under Mobius transformations, but as we will see one can define well-behaved 
operators with correlation functons that transform simply. 

The potential has a discrete shift symmetry N — > N + 1. The natural operators to 
consider have definite charge under the shift symmetry. The simplest such operators are 
exponentials, 

V p {z) = e i/3JV(2) . (3.4) 

We will see that these operators have positive definite weight at least under Mobius transfor- 
mations, and their correlators are finite in the IR. Correlators of N itself can be determined 
by differentiating the correlators of exponentials. As we will see this gives logarithms, as one 
would expect if iV were a massless field with dimension zero. Since the correlation functions 
of such fields are not well-defined, this is another reason to consider exponentials. 



3.1 The 1-point function 

To compute the 1-point function (Vp(zj) one simply needs to insert it into the partition sum 



w*»=*-Mi:£n 



n=0 k= 



dri 



\Js n 



d 2 x^e 



i0N+(z) 



^ n \ n 



n r R 



v n=0 



dri 



\Js n 



d 2 Xke 



-iPN-(z)) 



(3.5) 

where the two terms correspond to the instantons and anti-instantons. The exponential 
operator has a simple product form, and the contributions from the instantons and anti- 
instantons are complex conjugates of each other, so 

2 



(Vp(zi)) 



7-1 



exp 1 7+ ^ J d 2 x [e il3 B(2r - \x - Zx\) + 0(-2r + \x- Zi\)] J 



Cancelling against Z 1 and collecting terms, this becomes 

" R dr 



(Vpizi)) = exp 



-7(1 — cos (5) 



5 r " 



d 2 x&(2r- \x-zx\) 



(3.6) 



(3.7) 



The integral over x gives the area A\(r) such that a disk of radius r covers Z\ if its center is 
contained in A\. This set of points is a disk of radius r centered at z. Then 



(^(21)) = exp 



—7(1 — cos P) 



dr 



Mr) 



(3.8) 
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The area is A 1 = nr 2 , so the integral is 




Finally, the one point function of the exponential operator is 

{V ^ z))= [j) =\0 otherwise J" (3 ' 10) 

This type of "conservation of charge" condition is familiar from Liouville theory and free 
scalar CFTs, but the periodicity in (3 (which is a consequence of the quantization of N) is 
novel. 

3.2 The 2-point function 

One can compute the 2-point function by the same techniques. A similar analysis to the one 
above gives 



(VpM)Vp 2 (z 2 )) = \(exp[t(3 1 N(z 1 ) + t(3 2 N(z 2 )])\ 2 = (3.11) 
exp {-7 jf $ J d 2 x (1 - cos {PMr - \x - Zl \) + (3 2 Q{r - \x - z 2 \)})} . (3.12) 

This is equivalent to 

(VfsMW^z,)) = exp | - 7 J ^ [(1 - COS/30A? + (1 - cos/3 2 )A° + (1 - cos(/3i + P 2 ))A° 12 ] J 

(3.13) 

where we have defined the "exclusive area" A\ as the area of the region such that disks 
centered in that region cover z\ but not z 2 . Equivalently, construct two disks of radius r, 
one centered at Z\ and the other centered at z 2 . Then A® is the area covered only by the 
disk centered at z\. 

Similarly, we define the integral over these areas as 

I°( Zl ,Z2) = - ^A° 1 (r,z 1 ,z 2 ) . (3.14) 
7r J r 6 

With this notation, the two point function is 

(VpMWM) = exp {-7T7 [(1 - cos A)/? + (1 - cos/3 2 )/ 2 ° + (1 - cos(/3 1 + (3 2 ))I° 12 \ } 

(3.15) 

The integrals have a simple interpretation as spacetime volumes in de Sitter space: 

is proportional to the spacetime volume available to nucleate bubbles which cover only the 

points (zi,Zj,z k ). 
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To evaluate the integrals, we need to write the exclusive areas in terms of the simpler 
inclusive areas. For the two-point function, we have 





= A x - 


A 


A\ 


= A 2 - 


A 


A 

A 12 


= A 12 





(3.16) 
(3.17) 
(3.18) 



The same equations hold for the integrals of the areas, 



J? = h-I 12 (3.19) 
J 2 ° = J2-/12 (3.20) 
A° 2 = In (3.21) 

We need to compute Ji 2 , the integral of the area covered by both disks. Some simple geometry 
yields 

A 12 (r,z u z 2 ) = 2r 2 (cos^d^l^r) - (\z 12 \/2r)y/l - (\z 12 \/2r)^ 6(2r - \z 12 \), (3.22) 
where \z± 2 \ = \zi — z 2 \ is the distance between the centers of the disks. Integrating this gives 

h2{r,z u z 2 ) = - f R A 12 dr/r 3 = ln(R/\z 12 \) - 1/2 + 0(1/ R). (3.23) 
Plugging this in gives 

(3.24) 
(3.25) 



jO _ jO 
J l — J 2 


= In 


\zi 2 \ 
5 


1 

+ 2 


1° 
1 12 




R 


1 


= ln 


\ Z 12\ 


~ 2 



Let us redefine the UV cutoff 5 to eliminate the annoying constant factor so that 

It = I° 2 = In M (3.26) 

Note that I® and I 2 are infrared finite, while I u diverges as R — ► 00. This corresponds 
to an infinite expected number of disks covering both points 1 and 2. Therefore, for fii 7^ f3 2 , 



the exponent of Eq. 3.15 goes to —00 as the infrared cutoff R is taken to infinity. So the 



two-point function Eq. |3.11| is zero due to IR divergences. To cancel this divergence it is 
necessary and sufficient to require that the coefficient of the double overlap region If 2 is zero; 
in other words one needs cos(/9i + (3 2 ) = 1, or (3\ + (3 2 = 2im {n an integer). 

This condition is a kind of charge conservation condition: under the shift symmetry 
N — > N + 1, the operator Vp transforms as Vp — > exp(if3)Vp. So the correlators are nonzero 
only when they are invariant under the shift symmetry. Because iV takes integer values, 
the operator e %l3N ^ is equivalent to the operator e t W+ 2n ) N ( z ) ; so it is natural that the charge 
conservation condition is defined mod 27r. We will see in a moment that the dimensions of 
operators are also periodic functions of (3. 
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Enforcing this condition we get 



/ X \ 7^(1 - COsft) / r n 7T7(1 -cos A) 

<Vfc(*W*)> = (— ) (— ) (3.27) 

when ft + ft = 27rn (else the correlator is zero). Defining A(ft = A(/3) = §7(1 — cos/3), 
this has the form of a two-point function for a conformal operator of dimension (A, A). 



3.2.1 Correlation functions of N 

At this point, we pause for a moment in our analysis of exponential operators to consider 
quantities which may seem more basic: correlators of the field N itself. Because N — > —N 
is a symmetry, the 1-point function vanishes, (N(z)) = 0. The 2-point function can be 
obtained by differentiating the 2-point function of exponentials, 

(N(zi)N(z 2 )) = -——ty^Vfrfa)) (3.28) 



evaluated at Pi — Pi — 0. A convenient form of the correlator to differentiate is (3.15). 
Differentiating and setting p 1 = p 2 = gives 

(N{zi)N(z 2 )) = tt^ 2 = 7T7 (in A- - (3.29) 

This is the correlation function of a dimension zero field. To put it in a more standard form, 
the additive factor of 1/2 could be absorbed into the infrared cutoff and the prefactor itj 
could be absorbed into a field redefinition of N. The presence of the infrared divergence 
means that these correlators are not really well-defined. This is not surprising: since the 
theory has a symmetry N —+ N + 1, N is not a physical quantity. However, the exponentials 
we have been considering are physical and have well-defined correlators. Exactly the same 
issues arise for a free massless scalar in two dimensions. 



3.3 The 3-point function 

To compute the 3-point function of exponentials (V l 3 1 (zi)V^ 2 (z2)V^ 3 (z 3 )) we will need to 
evaluate the integrals of the overlap regions of three disks of equal size, each centered on a 
point Zi where the operators is inserted. The calculation proceeds along the same lines as 
for the 2-point function; omitting some details one obtains 

(VfkiziWkWVfriza)) = exp {-vr 7 [(1 - cos ft)!? + (1 - cos/3 2 )! 2 ° + (1 - cos/3 3 )! 3 ° + 
(1 _ C0S ( /3l + ft)) j0 2 + ( i _ cos(A + ^)) j0 3 + ( i _ cos(/?2 + ^ 3 )) 7 3 + (1 _ cos(/?1 + p 2 + ^ 3 ))/0 23 ] } 

Since the integral !? 2 3 is again logarithmically divergent at large r, a "charge cancellation" 
condition is required to cancel the IR divergence that would otherwise send the correlator 
to zero. Requiring the coefficient of this term be zero means ft + ft + ft = 2-7m. 
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When this condition is satisfied, the correlator simplifies to 



<Vk(*i)Vk(* 2 )Vfe(* 3 )> = 
exp {-7T7 [(1 - 008/30(7° + J 2 ° 3 ) + (1 - cos/3 2 )(/ 2 ° + i? 3 ) + (1 - cos/3 3 )(/ 3 ° + J? 2 )]} 

The exclusive area integrals are given by, for example, 

i? = II - /l2 - /l3 + Il23 (3.30) 
-^23 = -^23 — ^123 (3.31) 

The formula for the triple overlap /123 is somewhat complicated, but it cancels in the 3-point 
function, because the integrals appear in combinations such as 

A° + / 2 ° 3 = h~ I12 ~ lis + I23 (3-32) 
Therefore the 3-point function simplifies to 

exp {— 7T7 [(1 — cos — 1 12 — I13 + ^23)]} x (cyclic permutations) 
In terms of the weights A, = §7(1 — cosfli) the 3-point function can be written 

(Vp 1 (z 1 )Vp 2 (z 2 )Vp 3 (z 3 )) = |exp {-(Ai + A 2 - A 3 )(/ 1 - ii 2 )}| 2 x (cyclic permutations) 

(3.33) 

The combination (ii — Ji 2 ) is exactly the same as in the 2-point function, so 

/ X \ A1+A2-A3 / X \ A1+A3-A2 / X \ A2+A3-A1 

iVMVMVM) - (— ) (— ) (— ) 

(3.34) 

This is the 3-point function required by conformal invariance for three operators of weights 
Aj, Ai, with Aj = Aj. It is worth noting that scale invariance alone is not enough to fix 
this form — scale invariance requires only that the total scaling dimension of any term on the 
right-hand side be consistent with the total scaling dimension of the fields in the correlator, 
but not this particular structure. However global conformal invariance does require this form 
(in any number of dimensions), because there are no conformal invariants that can be built 
from less than 4 points. 



Starting from the Mobius invariance of the distribution Eq. 2.7 one could presumably 
prove that well-behaved correlators must be of this form. The statement is non-trivial 
because of the issue of IR divergences; correlators that depend on the IR regulator will not 
in general respect this form. 
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4 Four-point function 



In this section we find the 4-point function and analyze its properties. By now the procedure 
is familiar. The 4-point function is 

(V^V^V^V^z,)) = exp{-7r 7 E(l - cos &)/° + J> " cos (A + Pj)K+ 

i i<j 

(1 - C0S(A + ft + k ))I° k + (1 - COS Pi) A°234]} 

i<j<k i 

(4.1) 

The charge cancellation condition works as before: /1234 is the only infrared divergent quan- 
tity so the correlator is zero unless 

X> = 27m (4.2) 

i 

Using the charge cancellation condition, the 4-point function is 

exp j-A^/f + J 2 ° 34 ) - A 2 (/ 2 ° + 7 3 ° 41 ) - A 3 (/ 3 + J 4 ° 12 ) - A 4 (/ 4 ° + I° 123 ) - £ A y J§ J 

(4.3) 

where we have defined 

Ay = | 7 (l-cos(A + &)) (4-4) 
Note that due to the charge cancellation condition A12 = A34. 

In the case of four points, the exclusive area integrals are given by 



J? 


— h ~ I\2 — hs ~ I\& + A23 + I\2A + A34 — -^1234 


(4.5) 




= 7i2 — 7i23 — 7i24 + A234 


(4.6) 


7° 

i 123 


= 7/l23 ~~ A234 


(4.7) 


J 1234 


= -^1234 


(4.8) 



Using these relations and massaging the expression, the 4-point function can be rewritten 

(VpM)Vp 2 (z 2 )V P3 (z 3 )V P4 (z 4 )) = 



exp|-^(A i + A j -A ij )(7 1 -7, J )- (j2 A *-\j2 A v) (j2 ^ - 2I 1234 - 2I ± ^ 

I i<j \ i i<j / \i<j<k / 



(4.9) 

The first sum has the form of a product over six 2-point functions, while the second term 
contains a nontrivial function of the positions of the points. To be explicit, 

(V P M)Vp 2 (z2)V P3 (z 3 )V P M)) = 



U(j-) 3 'exp]- (E A *-^E A ^) ( E A,"* " 2/1234 - 2/i 

i<3 ^ lJ ' L \ i i<j J \i<j<k 
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(4.10) 



The interesting functional dependence on the location of the points is all contained by the 
function f(z u z 2 , z 3 , z 4 ) = Y,i<j<k Af* ~ 2/ i234 - 2J X . 

It is possible to compute the 4-point function in full generality and show that it is 
conformally invariant. However, since the distribution of bubbles is conformally invariant, 
the 4-point function is guaranteed to be conformally invariant unless infrared divergences 
arise. Therefore, we will assume conformal invariance and compute the 4-point function with 
the four points at 

zi = z (4.11) 
z 2 = (4.12) 
z 3 = 1 (4.13) 
z A = oo (4.14) 

With this assumption, the four-point function is an infinite constant times a nontrivial 
function of z, 



(^(^^(0)^3(1)^(00)) = C 



-(A 1+ A 2 -A 12 ) (1 _ z) -(A 1+ A3-A 13 ) exp 




It now remains to evaluate the integral J123. 



4.1 Evaluation of the triple overlap integral 

We need to evaluate 




(4.16) 



where j4 123 is the area contained within the triple overlap of three disks of radius r centered 
at the points zj, z 2 , and z 3 . 

To evaluate this we will need a formula for the area of triple overlap of three circles. In 
some cases this reduces to a double overlap, but in situations where the triple overlap is of 
triangular type (e.g. a region bounded by the arcs of three distinct circles) the area is [22J: 

^123 = \ [A 12 + A 13 + A 23 - vrr 2 ) + A T , (4.17) 

where At is the area of the triangle with vertices at the three points. We will do the 
computation for triangles for which the triple overlap area is always given by this formula 
for any disk size r, which amounts to assuming that the triangle is sufficiently close to 
equilateral, with none of the angles exceeding 90°. However, our final formula will be valid 
for any arrangement of three points. 

The triple overlap integral begins to be nonzero at the smallest value of r such that a 
disk can cover all three points. This value is called the circumradius R c . By subtracting the 
area inside a wedge of the circle from a triagular area (see figure |6| , we find 

i,4 2 3 = r 2 (n/2-9) - r 2 sin# cos 9 (4.18) 
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Figure 6: It is convenient to change variables in the integration, using 9 shown in the figure 
instead of r as the integration variable. 



where 9 is the angle shown in the figure. The integral of A 23 is infrared divergent, but the 
quantity A 23 — nr 2 is infrared safe, 

-(A 23 - vrr 2 ) = -{9 + sinflcos#)r 2 (4.19) 



so we focus on it. We need to integrate the triple overlap from the lower limit R c where it 
is first nonzero. So the A 23 piece of the integral is given by 



R 



dr 



{A 23 - vrr 2 ) 



Rc 



dr 



[9 + sin 9 cos 9} 



(4.20) 



and we can now freely take R — > 00 because the integral is finite in the infrared. The angle 
9 is related to r by 

sin#=— (4.21) 
2r v ; 

Also, as shown in figure [7j the lower limit r = R c corresponds to an upper limit on 9, 9 — 9%, 

where 9\ is the angle of the triangle with its vertex at point 1. The upper limit r = 00 

corresponds to 9 = 0. 

Making the change of variables, the integral is 



dr 



(A 23 - vrr 2 ) 



Rc 



d9 [9 cot 9 + cos 2 9] 



We can now write the full integral. 



Rc 



dr 



-A 



12:-! 



R 



dr 



((A 12 - ixr 2 ) + (A 13 - vrr 2 ) + (A 23 - Trr 2 )) + Trr 2 + A? 



(4.22) 



(4.23) 
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Figure 7: The smallest radius disk which can cover all three points is defined by r = R c . 
The figure shows that the lower limit if integration corresponds to 9 = 9\. 



Now the only infrared divergent term is the trivial 7rr 2 ; for the rest of the terms we take 
R — > oo to get 

/ —Ava = 7r tog — + {6 cot 6 + cos 2 9)d9 (4.24) 

JRc r K c ^ti c i JO 

where the sum is over the three angles of the triangle. 
The integral gives 

I \9 cot 9 + cos 2 9)d9 = ^ln(sin^) + 9 t {\ + ln2) + jsin(2^) + \% [U 2 (e 2l6 >)] (4.25) 
Jo 2 4 2 

where Li 2 is the dilogarithm function, 

oo n 

Li 2 (z) = J2 Z ~2 ( 4 - 26 ) 



7r 

n=l 



We would like to write everything in terms of the angles and lengths of the sides. By 
examining figure [7j we find 

2R C = A t = A t = A t (4-27) 

sin 9i sin 9 2 sin 6 , 3 



A T = R 2 C ^2 sin 9i cos 6i (4.28) 
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We can now rewrite the integral as 

tt f8R 3 sin 6»i sin 2 sin 6> 3 \ 1 ^ . 
3 V d 12 d 13 d 23 J 2 ^ 

\ ( 4 -29) 

-J] Uln(sin0,) + 6i(- +ln2) + ^ sin(2^) + ^ [Li 2 (e 2 ^)] 



In combining the terms, some simplifications occur because ^2 9i — tt. A nice symmetric 
way to write the answer is 



TT 7T i? S 
7^123 = -- + "In 



2 3 KL2^13^23 



£ (i6 t - ^)ln(sin^) + Icj [Li 2 (e 2 ^)]) (4.30) 



The sum is over all three angles of the triangle formed by the three points, and each angle 
is defined in the conventional way so that < Q{ < tt. Although our derivation is only valid 
for triangles which are sufficiently close to equilateral, the answer written in this way is valid 
for any arrangement of the three points. 

The answer can also be written in the pleasing form 

7r/ 123 = -| + tt In R - 9 1 ln d 23 - 9 2 ln d 3l - 6 3 ln d 12 - l - ^ 3 [Li 2 (e 2 ^)] (4.31) 

i 

It is particularly simple when the three points are on a line. If point 2 is between points 1 
and 3, then 9i = 9 3 = and 9 2 = tt, so that 

J 123 = -1/2 + InR - \nd 31 (collinear) (4.32) 

because all of the dilogarithms vanish. 

Having written the answer in terms of the angles, we want to write it as a function of z 
for our special choice of points, 

Z\ — Z 

z 2 = 

(4.33) 

The angles should satisfy < 9^ < tt. For $$(z) > we have 

e« = ,4.34) 

e 2 ^ = - (4.35) 

z 

e 2 ^ = i— ^ (4.36) 



18 



The above is not valid for ^s(z) < 0, but it is clear that the correct answer should be invariant 
under z — > z, so we will just work it out for z in the upper half plane and determine the 
value in the lower half plane by symmetry. 

We rewrite the answer piece by piece in terms of z. To write the dilogarithm part of the 
answer in terms of z, it is convenient to rewrite the answer in terms of the Bloch-Wigner 
function D2. The relation is 

D 2 (z) = $t[U 2 (z)] + arg(l - z) In \z\ (4.37) 
so that Y) 2 {e %e } = SfL^e )]. The combination which appears in the answer is 

\ E ^- W = H D2 © +D2 (§^I) +D! (^)} = (438) 

where we have used an identity of the Bloch-Wigner function to get the last equality (see 
section 7.2). 

For the part not involving dilogarithms, we need to solve for 9i in terms of z. Again for 



z in the upper half plane, we can invert (4.36) to get 

%d 2 = Infj^i 

-iOz = ln (|f3^) ( 4 - 39 ) 

subject to the usual convention that the branch cut for the logarithm is taken to be on the 
negative real axis. 61 is not needed because in the answer it multiplies lnc^, which is zero. 
Plugging these in and simplifying, 

- 6 l In d 23 - 6 2 In d 3 i - 63 In d 12 = -3 [In z ln(l - z)} (4.40) 

So the answer as a function of z is 

7rI 123 = -| + 7rlni2-3[lnzln(l-;z)]-D 2 (z) for %z) > (4.41) 

Now our algebra has been done under the assumption that z is in the upper half plane, 
but the answer must be symmetric under z — > z. The answer as it stands is a constant 
term plus two functions which are odd under z — > z. Therefore, the correct answer has an 
additional factor of the sign of the imaginary part of z, 

Trim = -^ + 7r\nR- sgn(3(z)) {3 [In z ln(l - z)\ + D 2 (^)} (4.42) 



4.2 Properties of the 4-point function 

Plugging this in, we have the 4-point function 

(V Pl {z)V P MV^{l)V^)) = C \z- A ^ +A ^(l - z )-a 1 -a.+a u |2 x 



exp 



lJ2 Ai ~^J2 Ai n s ^( z )) {^t ln z ln ( x - «)] + D 2 } ( 4 - 43 ) 
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By conformal invariance, aside from simple prefactors which depend on the conformal weights 
of the fields, the 4-point function can depend only on the cross ratio. For our choice of points 
the cross ratio is just z: 

{z\ - z 2 )(z 4 - z 3 ) 



(z 3 - z 2 ) (z 4 - Zx) 



(AAA) 



The 4-point function in a conformal field theory should be crossing-symmetric — it should 
be invariant under interchanging the points where the operators are inserted. Interchanging 
the points corresponds to the following group of transformations on z: 

111 —z 

z -> 1 - - -> — ► - -> 1 - z -> — (4.45) 

z 1 — z z 1 — z 

The last three are odd permutations of the four points; the rest are even [25]. The Bloch- 
Wigner function D 2 (z) has the property that it changes sign under odd permutations and is 
invariant under even permutations j2S] . Since the sgn function changes sign only under odd 
permutations, sgn(Q ; (;2:))D2(z) is invariant under all permutations. 

The other nontrivial factor in the 4-point function, (sgn(^(z))S[ln2:ln(l — z)]), is slightly 
more complicated, but one can check that under permutations it transforms in the correct 
way to contribute the right factors. To take a nontrivial example, under z — > 1/z this 
function transforms as 

sgn(S(z" 1 ))3[lnz- 1 ln(l - z~ 1 )} = sgn(9f(z))9fpiizlix(l - z)) - Trln \z\. (4.46) 

One can use this to check that the 4-point function Eq. |4.43| satisfies 

(^(^(0)^3(1)^(00)) = ^^(^(1/^)^(0)^3(1)^(00)}, (4.47) 

which is the correct behavior [21] . 

The four-point function simplifies when all 4 points are on a line, or more generally a 



circle. Using the simple form (4.32) for the triple overlap integral, for z real and negative 
we have 



(V A (Z)V A (0)V A (1)V A (00)> = C Is-Ai-^+A^ _ ^-Ai-Aa+A^ x 

exp < 2 Aj - Ay ) In |1 - z\ > (z real and negative) . (4.48) 




Therefore for all 4 points on a circle, the 4-point function is simply a product of power laws; 
the powers which appear depend on the order of the points. 

However, there is a problem. Without the sgn factor, the functions in the exponent in 



Eq. |4.43| would be odd under z — > z. Furthermore, they are real analytic functions away 
from lines of discontinuity running along the real axis from —00 to and from 1 to +00. So 
before multiplying by the sign function, the exponent is a real analytic function in a finite 
region around z = 1/2. Therefore after including the sign the 4-point function is not real 
analytic as we drag z across the real axis near z = 1/2, even though z is separated from 
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the other points at and 1 by a finite distance. In a conventional field theory, correlation 
functions should be analytic except when two points approach each other. 

The nonanalyticity of the 4-point function indicates that the system we have defined does 
not correspond to a full-fledged conformal field theory. This is not too surprising, because 
we have made a number of approximations in treating the physics of bubble nucleation; one 
may still expect that our system can be obtained as the limit of a genuine conformal field 
theory. We will discuss this further in the conclusions. 



5 Central charge 

Conformal field theories on curved spaces have a conformal anomaly. Specifically, in 2D 
CFTs T a a = —j^R, where R is the scalar curvature of the 2D space. Since the trace of the 
stress tensor is related to the variation of the action with respect to the conformal factor in 
the metric, one can compute c by taking the derivative of In Z with respect to the log of the 
curvature. On a sphere of radius R, 

Z(R) = R c/3 Z 0l (5.1) 

where Zq is the partition function on a sphere of unit radius [23J. 

Given the results of the previous section, it is at best unlikely that the theory as we have 
defined it is a full conformal field theory. Nevertheless we will proceed, as we can easily 
compute Z on a sphere. The partition function is 



n=l 



where e is a cutoff on disk angular size and A is a (cosmological) constant added as a local 
counterterm to cancel the leading UV divergence from small e (as a constant multiplicative 
factor in Z it cancels out of all the correlators computed earlier). Computing the integral 
and restoring the dimensions gives 



In Z = 47T 7 (% + In - + In 2 - - - AR 2 + 0(e 2 )] 
\ e 2 e 6 / 



(5.3) 



Setting the counterterm A = 1/e 2 cancels the quadratic UV divergence, but the log is an 
anomaly that cannot be cancelled with any local counterterm. This is precisely what one 
expects for a 2D CFT on a sphere with central charge 

c = 12tt7 . (5.4) 

For c < 1, unitarity implies that c must take a discrete set of values (the minimal models). 
Therefore if this calculation is taken seriously it indicates that our model cannot be unitary 
at small 7. However it is worth mentioning that we have included neither perturbative 
fluctuations of the field nor of the geometry. One expects graviton fluctuations to contribute 
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a term of order (Mp/H) d ~ 2 , and so it is possible that the term we have computed is only 
one of several contributions to the central charge of a putative complete theory. At large 7, 
it is reasonable to expect that the instanton fluctuations are the dominant contribution, and 
indeed in that limit the theory becomes free. It would be interesting to compare this limit 
to the analogue in AdS space, where one would take Gn — > with the AdS radius held fixed 
so that c — > 00. 

5.1 CFT on a fractal 

In the string theory landscape there are so-called "terminal vacua" : minima with either zero 
or negative cosmo logical constant. A region which tunnels to one of these minima has at 
most a finite probability of nucleating any more bubbles before infinite time (for a zero CC 
bubble) or a big crunch (for a negative CC bubble). How best to deal with these regions 
is unclear, but the proposal of [13] is to excise them and attempt to define a CFT on the 
remaining space, perhaps including lower dimensional defect CFTs on the boundaries. 

A simple toy model for terminal vacua in our 2D CFT is to assume there is some rate 
7 to produce "dead" disks. Following the suggestion of [T3] then corresponds to computing 
correlators in the regions outside these dead disks. Since the region covered by zero disks of 
any type is a fractal set of measure zero, these correlators will "live" on a fractal set. 

One immediate problem is that when 7 is large this set is not just measure zero, but 
empty (recall that dp = d — 7^2^, and when dp < the set is empty). One can consider the 
case of 7 < 1 and continue, but it is clear that the resulting theory will have a very different 
structure than what we have considered so far. 

Since the set is by definition N(z) = one cannot compute correlators of the disk number 
operator in the way we have been proceeding. There are two obvious approaches one could 
take to this. One is to compute the probability that some set of points Zi in the full space 
are all in the set, or various conditional probabilities (such as the odds that if one point 
is in the set, the others are as well). Computing these probabilities is not difficult using 
the techniques we have already developed, and the results depend on the distances between 
the points. However proceeding in this manner we have not succeeded in defining a set of 
probabilities that are independent of the IR cutoff. The problem in a nutshell is that the 
number of disks which covers some but not all of a certain set of points is IR finite when 
integrated against the distribution, but on the other hand the number of disks that would 
cover all of the points is infrared divergent. 

Another approach is to consider more types of disks: a "dead" type which defines the 
fractal, and then one or more other "live" types. One could then try to compute correlators 
of the number operator for "live" disks within the set of points covered by zero "dead" disks. 
However this analysis requires taking into account interactions between the different types 
of bubbles, something which we will not consider in this note. 
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6 Generalization to arbitrary dimension 



The definition of the theory can be easily generalized to arbitrary dimensions. The partition 
function in d dimensions is 

Z = Z + Z_ (6.1) 

with 

z+ =^{^Ss^S ddx ) (6 - 2) 

We will show that the two and three point functions generalize in a simple way. To make 
the discussion here easier to follow, we include some formulas and discussion which overlap 
with the earlier sections of the paper. 

The general N-point function is 

n 

(jthNMjfhNM^i&NM} = exp(-2C d7 [J](l - cos + $^(1 - cos(ft + + 

i=l i<j 

^(l-cos(A + ^ + A))4 fc + ...]!6.3) 

i<j<k 

where now Zi is a point in rf-dimensional space and C d is the volume of a unit d-sphere. As 
before, define the integral 

^/^M- (6-4) 

The two-point function when the charge conservation condition is satisfied is 

( e tf"(*i) e -*/w(«)) = exp(-C d7 (l - cos/3)(7 1 ° + J 2 )) (6.5) 
Now I® = I® = h — I V2 as before. The 2-point function is then 

{e iPN le - iP N 2) = exp( _2C d7 (l - cos (3){h - J 12 )) (6.6) 
The integral is given by 

A - /i2 = 7T / -^Mr) - / -fZi^(r,d 12 ) (6.7) 



C d J s r*" 1 ^ > C d J dl2/2 r^ 



It is helpful to rewrite this as 



The first term can be integrated immediately using A\ (r) = Cd,r d to get 

h - In = In (I) - i- £ £ [*(r) - *,(r, «,„)] (6.9) 



23 



The remaining integral can be evaluated explicitly, but the crucial information can be ex- 
tracted more cheaply. The first term, the integral of Ai, diverges logarithmically at large 
r. However, the second term cancels this divergence, because the double overlap region A 12 
asymptotically has the same area as a single disk, 

A 12 -> C d r d as r -> oo . (6.10) 

Therefore, the integral is infrared finite, and we can take R — > oo, so that the integral is 



1 



dr 



C d J*i2 r^ 1 



[A(r)-A 12 (r,rf 12 )] 



(6.11) 



This formula depends only on the dimensionful quantity du, but dimensional analysis shows 
that the answer must be dimensionless. Therefore it is a constant independent of the distance 
di2- So finally 



h-I 



12 



In 



12 



d_ 

~5 



- #d 



(6.12) 



where is a dimension-dependent constant. 
Plugging this in, we have 



(e 



iPN{zi)-iPN{z2) 



-2C d7 (l-cos/3)# d 



F12 



2C d7 (l-cos/3) 



(6.13) 



The prefactor can be absorbed into a redefinition of the ultraviolet cutoff 8. The dimension 
of the vertex operator e tf3N is 



A = A 



C d 



7(1 — cos (3) 



(6.14) 



Having found the dimensions we can rewrite the general n-point function 



^ e i(3iN(zi) e il3 2 N(z2) e i/3 n N(z n )^ _ 



expj-f>/°- £ A yfc /° fc -... 

L i=l i<j i<j<k 



(6.15) 



Now to evaluate the 3-point function explicitly. In general, we have 

(e ^(*) e ^N(«) e ^(«) > = | exp |_ Ai/ _ A2 j0 _ As/ _ Ai2/ 2 _ Ai3/ 3 _ A23/ 3 _ Al23/ 23 ]}p 

(6.16) 

unless A123 = because is infrared divergent. A123 = 
londition is satisfied, Pi + P2 + P3 = 27m with n an integer. 



— ±40 

(6.16) 

The correlation function is zero unless A123 = because J° 2 3 i s infrared divergent. A123 = 
when the charge conservation condition is satisfied, Pi + P2 + P3 = 27rn with n an integer. 
Using the charge conservation condition, we have relations like A 12 = A 3 . Also, recall that 
the exclusive volumes are given by 

— Il3 ~ h3 + 1 123 



I°1 = h-hl 
Il2 = 1 12 ~ h23 



(6.17) 
(6.18) 
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so the 3-point function becomes 

|exp {-Ai(/i - /12 - /is + hs) ~ A 2 (/ 2 - In - hs + hs) - A 3 (/ 3 - I32 - hs + /i 2 )])}| 2 

(6.19) 

Note that the triple overlap region 7 123 does not appear. Collecting terms, we have 

|exp {-(Ax + A 2 - A 3 )(/i - /12) - (Ax + A 3 - A 2 )(h - J 13 ) - (A 2 + A 3 - A 1 )(/ 2 - J 23 )}| 2 

(6.20) 

where we have used I\ = J 2 = J 3 . But this factorizes into 2-point functions! This is precisely 
the form the 3-point function must take due to conformal invariance. Explicitly, it is 

r2 r2 r2 

^ e iy9iJV(zi) e i/3 2 N(z 2 ) e ip 3 N(z 3 )j _ ^ ^A!+A 2 -A 3 ^ ^Ai+A 3 -A 2 ^ ^A 2 +A 3 -Ai 2 1) 

^12^12 ^13^13 ^23^23 

where we have absorbed the same constant factor into the definition 8 as in the 2-point 
function. 
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Now for the 4-point function. Performing a similar analysis as for the 3-point function, 
we find 

2 

(e iA Nl e i/32N2 e iA7V3 e iPiNi ) 

i<j 

|exp {-Ca(/i 23 + /124 + /i34 + ^234 - 2/i 234 - 2h)}\ 2 (6.22) 

with the definition 

Ca = Ai + A 2 + A 3 + A 4 - Ai 2 - Ai 3 - A u ■ (6.23) 
The first factor has the form of factorized two point functions, so we can rewrite this as 

( e i(hm e ifhN ae ifkN 3e ifUN^ = jj / ° \ x (6 24) 

iKj Vl%l / 

|exp {-C A (/ 123 + 7i24 + /i34 + / 23 4 - 2/ 1234 - 2/i)}! 2 (6.25) 

So the first part of the 4-point function consists of simple power laws. The last line has all of 
the interesting information in it, and involves the triple and quadruple overlaps. This term 
is IR finite on its own, but it does depend on the UV cutoff through Ji; this dependence is 
trivial. 

Computing the 4-point function explicitly is a nontrivial task which we have only accom- 
plished in d = 2 so far. 



6.1 Free field limit 

Our correlators have a free-field limit when the tunneling rate becomes large. In taking 
7 — > oo, the expectation value of N will become very large. Operators which are well- 
defined in this limit should have fa — > 0. More preciesly, defining 

p = a J— (6.26) 

V 7T7 
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we want to take the limit 7 —>■ 00 with a fixed. 



To see this that the correlators factorize in this limit, start from the formula for an 



n-point correlator in general dimension, Eq. 6.15 



l e i(3iN{z 1 ) e ip 2 N{z2)^ Ji ip n N{zn)\ 



exp \ - A,i? - £ Ay/j - A ^4 k - •• 

i<J i<j<k 



i=l 



The exclusive integrated areas satisfy the obvious generalizations of Eqs. 



(6.27) 



4.5 



n ~ k (—1)' 

A°ii 2 ...ifc = ~~ /j liin-ikh-h (6.28) 

'=0 ' 3l,32,-,3l 

where the j indices are summed from 1 to n and the factor of Z! corrects for overcounting. 
In this expression we have defined hjk... = if any i,j,k... are equal. Then the exponent in 
Eq. |6.27| can be written 



fc! Z^ l\ Z^ J *X*2-»W1-JI 

k=l ' 1=0 ' jl,32,-,3l 

where we have also defined all areas with m indices Ii 1 i 2 ...i m = if m > n. 



(6.29) 



To continue, we would like to evaluate the coefficient of the areas Iixi2...iiji...j m -i with vn 



indices, which is: 



m— 1 



y a- ■ V I 



Z!(m-Z) 

z=o v ' 

In the limit j3i — ► 0, we have 



12— «m-iJl— 



Jl ,32,— ,3l 



Hi. ..Ik 



tt 7 (1 - cos^ + ... + AJ) = K + ... + «ij 2 + C(a 4 /7)- 



Since i^fe... is completely symmetric, 



A ii...i m -i^iii2:.i m -iji.:ji 



((m — Z)a^ + (m — Z)(m — I — 1)0^0^) 1^ 



I2...i m -i3\...3\ 



(6.30) 

(6.31) 
(6.32) 



Therefore Eq. |6.30| is equal to 



y , . y 1 



il,...,i m =l 



Z=0 



«7 



+ 



Z! V(m-Z-l)! (m-Z-2 



' ii,...,i m =l 



(6.33) 



So only the 1- and 2-disk overlap areas contribute! Therefore the exponent in Eq. |6.27 
becomes simply 

— afli — atiajlij. (6.34) 

Recalling that U = \n(R/S) and = ln(i?/|^-|)(l— dij), and using £\ a l = ( a i a j ~ ^^(l 
we get 

2 



JfaNizi) e ip 2 N(z2) e i)3 n N(z n ) 



exp j - ^ a i a 3 y + ln |~| (! ~ <%)J 



26 



(6.35) 



As usual this is IR divergent unless we impose conservation of charge, which here is simply 
Tlii a i = T^ij a i a j = 0- Then, recalling that (3N{z) = y/2a(j>(z) 

(6.36) 

This is the general form of a correlation function of vertex operators e l ^ a ^ z ^ of a free 
massless field in D = 2 (see e.g. [2T], p. 296). In higher dimensions a non- interacting 
massless scalar with non-canonical kinetic term (pO d ' 2 (p could produce such correlators. 



( / e iv / 2ai</'(2i) e iv / 2a2</>(z2) e i\/2a n <t>(z n )\^ _ J~J 



7 Conclusions 



Starting from eternal de Sitter space we have successfully defined a model with correlation 
functions that are conformally covariant and transform with positive weight. However be- 
cause of the lack of analyticity in the 4-point function, the model does not appear to be a 
healthy conformal field theory (except perhaps in the non-interacting 7 — > 00 limit). 

One possibility is that the theory we have defined here is the limit of some good CFT 
in which certain effects have been ignored. Such limits can result in non-analyticities in 
correlation functions (for example, one can get logs from power laws by expanding around a 
limit where the dimensions of some operators go to zero). Adding weight to this possibility 
is that in defining the simplest possible non-trivial model, we indeed have ignored many 
potentially important effects: 



• We have ignored perturbative fluctuations of the field around its minima, and included 
only the instantons. 

• We have ignored perturbative corrections to the instantons themselves, which for ex- 
ample include fluctuations away from spherical shape. 

• We have ignored interactions between the instantons other than their collisions, and 
we have treated collisions and overlaps in a simplistic manner. 

• We have used a semi-classical approximation that ignores quantum interference be- 
tween different configurations in the ensemble of bubbles. 

• We have ignored gravitational fluctuations in the bulk, as well as fluctuations in the 
geometry of the boundary. 

• The bulk theory we considered is a single scalar field in de Sitter space. Even coupled 
to gravity, such a model in anti-de Sitter space probably does not define a consistent 
CFT — one presumably needs the infinite number of modes of string theory, or at least 
one expects some very restrictive conditions on the bulk degrees of freedom necessary 
to define a good dual. Perhaps similar restrictions apply here. 



27 



From this point of view the model described here is only a potentially interesting first 
step on the path to a full theory of inflation or de Sitter space. Nevertheless, we feel the 
basic structure may be correct, and indeed many of the effects mentioned above could be 
included as perturbations around our limit. In addition, the model may have interested 
applications to condensed matter systems or as an example of a new class of conformally 
invariant theories. We hope to investigate some of these issues further in the future. 
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